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ABSTRACT 


A bow wave has been observed in the solar wind upstream of each of 
the six planets visited by spacecraft. The observed properties of these 
bow waves and the associated plasma flows are outlined, and those 
features identified that can be described by a continuum 

magnetohydrodynamic flow theory as opposed to a more detailed 

multicomponent particle and field plasma theory. The primary objectives 
of this paper are to provide an account of the fundamental concepts and 
, current status of the magnetohydrodynamic and gas dynamic theories for 
solar wind flow past planetary bodies. This includes a critical 
examination of (a) the fundamental assumptions of the theories, (b) the 
various simplifying approximations introduced to obtain tractable 
mathematical problems, (c) the limitations they impose on the results, 
and (d) the relationship between the results of the simpler gas 
dynamic-frozen field theory and the more accurate but less completely 
worked out magnetohydrodynamic theory. 

Representative results of the various theories are presented and 
compared. A number of deficiencies, ambiguities, and suggestions for 
improvements are discussed, and several significant extensions of the 
theory required to provide comparable results for all planets, their 
satellites, and comets are noted. The paper concludes with some remarks 
about anticipated trends in the future development and application of 
the theory and a description of a number of applications and extensions 
currently under development. 
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INTRODUCTION 

♦ 

Direct observations have established the presence of a bow shock 
wave upstream of each of the Innermost six planets of the Solar 
System. These bow waves form In the solar wind plasma* which flows past 
the planets at supersonic speeds carrying solar plasma and magnetic 
fields away from the Sun, and define with certain exceptions the 
upstream boundary of the region influenced by the presence of the 
planetary body. It is the purpose of this paper to review the 
quantitative magnetohydrodynamic and gasdynamic theories for calculating 
the properties of these bow wayes and their associated flow and magnetic 
fields. 

The principal features of these flows are as follows. Hot plasma 
of the sclar corona accelerates to supersonic speeds within a few solar 
radii of the Sun and flows outward through the Solar System. Because of 
.its large scale and high electrical conductivity, this plasma carries 
the embedded solar magnetic field with it, distorting it as it proceeds, 
and generally forming a spiral magnetic field configuration as a result 
of the comDined action of the outward flow of the plasma and the 
rotation of the Sun. At the orbit of Earth, the free-stream velocity v^ 
may be less than 300 km s“ l or greater thanlbOOkm s" 1 , the proton 
number density Np may occasionally be less than 1 cm -3 or as high as 
100 cm -3 , the proton temperature T p may be less than 50,000 K or as 

' oo 

much as 500,000 K, and the magnetic field B ra may be less than 1 nT or as 
much as 50 nT. The ion composition is primarily protons, i.e., ionized 
hydrogen atoms; with from 5 to 20 percent of helium and other ions. 
There is also an equal number of electrons, and their temperatures are 
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usually several times greater than that of the Ions. The solar wind 
flow approaches the Earth within a few degrees of the direction to the 
Sun, with the aberration due to the planet's orbital motion providing a 
substantial part of the deviation. The interplanetary magnetic field 
direction is particularly variable, with the greatest predominance being 
approximately in the plane of the ecliptic with an angle to the 
Sun-Earth line of approximately 45°, and directed about equally either 
toward or away from the Sun. 

Figure 1 illustrates the nature of the resulting flow about the 
Earth and other major obstacles in the Solar System' in barest outline. 
The Earth's magnetic field, being relatively strong, dominates the 
plasma in an extended region surrounding the Earth to form the 
magnetosphere. This region is bounded by the magnetopause, a relatively 
thin layer In which electrical currents flow in such a way as to bound 
the magnetospheric magnetic field on one side and the rapidly flowing 
solar wind plasma and magnetic field on the other. The solar wind 
plasma experiences considerable difficulty in trying to cross the 
magnetopause, and with some exceptions is excluded from the 
magnetosphere and forced to flow around the magnetopause much as though 
it were an impenetrable obstacle. The solar wind flow approaches the 
planet at speeds of the order of 5 to 10 times the fastest wave 
propagation speeds involving the ion motions, and decelerates abruptly 
through a bow shock wave much as air does as it approaches a blunt-nosed 
obstacle at supersonic speeds. Near the nose of this bow wave, the flow 
is slowed to subsonic speeds, and turned to begin its flow around the 
magnetopause. It accelerates as it flows, and attains supersonic speeds 
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as it proceeds toward the flanks about as sketched. For representative 

conditions in the solar wind, the geocentric distance to the 

» 

magnetopause nose is about 10 Earth radii, and that to the bow wave nose 
is about 3.5 Earth radii more. 

It is apparent that many features of this flow resemble those of 
familiar supersonic gasdynamic flow past a blunt-nosed obstacle such as 
a round-nosed bullet or a re-entry space vehicle. But gasdynamics alone 
is not sufficient to describe other important features of the flow, such 
as the formation of the magnetospheric obstacle itself, that depend on 
the effects of a magnetic field in a highly conducting plasma. The 

simplest theory that is capable of representing the main features of 

these flows is thus magnetohydrodynamics of a compressible highly 

conducting plasma. 

Although it is clearly established that magnetohydrodynamics can 
provide a remarkably good prediction of many features of these flows, it 
should be noted that this theory is not only complicated and far from 
completely worked out, but also fundamentally incapable of describing a 
number of other important features of the interaction. These include 
waves and energetic particles that sometimes propagate upstream from the 
bow wave along the interplanetary magnetic field, the differing nature 
of the magnetic field jumps across the bow wave depending on the angle 
between the incident interplanetary magnetic field and the shock normal, 
and the nonMaxwellian character of the proton velocity distributions. 
Analysis and prediction of such properties of these flows requires use 
of multi -component plasma theory as opposed to the single-fluid 
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magnetohydrodynamic theory under discussion in this paper; and is the 
subject of a number of other papers at this conference.' 

The Earth may serve as a prototype for discussing planetary bow 

shock waves and their associated flow fields, but there are also notable 

differences amongst the planets as outlined in Figure I. These stem 

from differences in the planetary magnetic fields and atmosphere, and 

from changes in the solar wind with distance from the Sun. As for the 

solar wind, the velocity remains virtually the same from the orbit of 

Mercury to beyond Saturn, the temperature, density, and- magnetic field 

* 

all decrease with distance from the Sun, and the sonic and Alfven Mach 
numbers of the flow incident upon the planets remain much greater than ' 
one; but the planetary properties differ quite markedly. Compared to 
the Earth, which is essentially a magnetic obstacle in the solar wind, 
Mercury has such a weak magnetic field and tenuous atmosphere that the 
distance r 0 to the nose of the magnetopause is only about 1.5 times the 
radius rp of the planet; whereas Jupiter and Saturn have such large 
magnetic fields that the ratio r 0 /rp has values of about 70 and 22 
respectively. The Earth's magnetopause is approximately axisymmetric 
about a line through the Earth's center and parallel to the free-stream 
direction; whereas the magnetopause of Jupiter, and perhaps of Saturn, 
is flattened significantly toward the equatorial plane of the planet, 
probably because of its rapid rotation and large size. Venus, which 
formerly seemed so much like a sister planet of Earth but is now known 
to be very different in many respects, has virtually no detectable 
magnetic field, but does have a sufficiently dense ionosphere to 
withstand the solar wind and form an obstacle having a shape as 
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Illustrated. In recognition that the boundary of the flowing solar wind 
bounds the ionosphere and not the magnetic field, it is called the 
ionopause. Its nose is only a few hundred kilometers above the planet's 
surface, and the nose of the bow wave is about a third of a planetary 
radius above that. The case of Mars remains uncertain because of the 
scarcity of observational evidence and the possibility that the location 
of the magneto ionopause can be accounted for by either a magnetic field 
or an ionospheric interaction. It is probably impossible to reach a 
definitive conclusion with the available magnetometer data, but other 
evidence indicates that there is only a small range of possibilities for 
the magnetic field that is compatible with the interpretation that Mars 
acts like a magnetic obstacle. If the magnetic field were substan- 
tially larger than the proposed values, the magnetosphere would be notice- 
ably larger. If the magnetic field were substantially smaller than sug- 
gested, it would have negligible effect on the location of the magneto- 
ionopause and the bow wave, and the obstacle would be clearly ionospheric. 


While not a planetary body, it is of interest to comment on the 
case of the Moon, whose environment has been explored extensively with 
plasma probes and magnetometers. Its interaction with the solar wind is 
quite different from that of the planets. Having neither an atmosphere 
no r significant magnetic field, the Moon is unable to stop the approach 
of the solar wind until it is absorbed at the lunar surface. The solar 
wind thus flows unimpeded into the sunward surface of the Moon, and 
leaves a cavity in the solar wind extending downstream from the night 
side. Instead of forming a bow shock wave upstream of the Moon, a 
trailing expansion Mach wave fan extends from the vicinity of the lunar 
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terminator as illustrated through which the solar wind plasma is first 

turned toward the axis to begin to fill in the cavity. 

* 

Still another class of related objects in the Solar System for 
which similar questions are being addressed are comets. At this 

writing, no measurements haye been made in situ at comets, but several 
spacecraft are already launched on their courses to rendezvous with 
comets in the near future. There is considerable uncertainty about the 
precise nature of the interaction, but the general consensus is that the 
situation is as illustrated in Figure 1. The comet nucleus -is a source 
of gas and particles that melt or evaporate off as the comet travels 
along the more sunward part of its trajectory. These particles, emitted 
as neutrals, become ionized as they are exposed to sunlight and the 
passing solar wind particles, and in effect appear to suddenly 

materialize in the ionized gas flow the moment they become ionized. 
This simulates addition (or subtraction) of mass, momentum, and energy 
in the flow; thereby requiring further fundamental considerations in 
the formulation of the theoretical model. If the ionization occurs 
sufficiently rapidly and near the cometary nucleus, the flow field may 
resemble that of the planets with an ionopause that separates most of 
the cometary material from the solar wind. If the ionization extends 
over a substantially larger region, the approaching solar wind flow may 
be altered considerably, even to the extent that there might be no bow 
wave whatsoever as the mass-momentum-energy addition slows and heats the 
oncoming plasma so as to gradually change the supersonic flow to 
subsonic. 
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Finally, there are the natural satellites or moons of the other 

planets. Not a great deal Is known about the flows about these objects, 

» 

although It Is evident that a number of statements can be made on the 
basis of their environment, sometimes in the solar wind upstream of the 
planetary bow wave, and sometimes in the planetary magnetoionosphere or 
the intervening magnetosheath flow*, and also on their own atmospheric 
properties. From the standpoint of magnetohydrodynamic and gasdynamic 
theories of planetary bow waves, the differences among the various 
planetary flows are to seme degree secondary as they are concerned more 
with the specification of appropriate boundary conditions than with the 
goyerning physical laws and equations. The remainder of this paper will 
concentrate on providing an account of the theoretical formulation and 
results of a number of mathematical models for these flows, and on 
outlining the various approximations that have been introduced to 
surmount the almost intractable computational problems posed by the 
magnetohydrodynamic model. 


MAGNETOHYDRODYNAMIC THEORY FOR PLANETARY BOW WAVE FLOWS 

The objective for theoretical studies of planetary bow wave flows 
i« to predict the bow wave location and the properties in the 
magnetoionosheath region of the associated flow field given the 
properties of the incident solar wind plasma and the planetary body. 
The simplest theoretical model that can provide the salient features of 
these flows is that provided by magnetohydrodynamics of a compressible 
highly conducting perfect gas. This theory, which has been gradually 
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developed over the last forty or fifty years represents a combination of 
single fluid, continuum gasdynamics, and Faraday electricity and 
magnetism, i.e., Maxwell’s displacement current is omitted on the basis 
of its smallness ir appropriate applications, This theory is most 
frequently presented in its perfect dissipationless form, but it is 
advantageous for the present discussion to include the dissipative terms 
representing effects of viscosity, electrical resistance, and thermal 
conduction. These terms will be written in the usual way (see e.g., 
Landau and Lifshitz, 1960, or Jeffrey, 1966) with scalar coefficients 
for first and second viscosity y and ?, electrical conductivity cr , and 
thermal conductivity k; but it should be recognized that the transport 
properties of a nearly collisionless plasma with embedded magnetic field - 
may be highly anisotropic with very different values parallel and 
perpendicular to the magnetic field. These considerations lead to the 
following set of partial differential equations for the pressure P, 
density p, temperature T, velocity v, magnetic field B, and entropy s 

*v 

involving the universal gas constant R = 8.31 x 10 7 ergs gm" l K" x , the 
mean molecular mass ffl of the gas in atomic units, and the gravitational 
potential $ and constant G = 6,67 x 10" 8 dyne cm 2 gm" 2 ; and in which 
the independent variables are the time t and space coordinates x of an 
inertial system. 


Perfect gas 

P = pRT/m 


( 1 ) 



Conservation laws 
Mass 


If + 3x£ (pv k J " 0 


( 2 ) 


Momentum 


3tt. 


JL( PVj ) 0 

3t ' i' 3X k 
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where 


B l B k , B 2 


«1*k 
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3v t 


9 i s " 3xJ * 5 a ik = y ^ax k + ax. ’ T 6 ik ax £ +?(S ik ax 5 
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)(5) 


Energy 


JL 'p y 
at 


(£ 75 — +06 + 0 $+ gjp) + div q s 0 


(6) 


where 


q = pv(-7j- + e + ~ + $ ) + B x (v x B) 


m; B x curl B-v^a'-kVI 

1 DTT 0 - — ** ** « ** 


(7) 


or the Heat Transfer Equation 


av 


pT If ■ «i' k ^ + div(k ? T) + raw (cu - rl ? )2 = div(k ? T) 

+ f + Y^> 2 + ? < div x> 2+ TfW< su - rl B)2 w 
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and the Faraday (preMaxwell) Electricity and Magnetism Equations 
33 r 2 

TO ■ c “rl (V X B) + Isy V*B , dlv B - 0 (9) 

In these equations tt^ and are the momemtum flux and viscous 
stress tensors, g^ is the acceleration of gravity, 6^ is the Kronecker 
delta function, q is the energy flux vector, and e * cjr is the internal 
energy. Other relations not explicitly used in the above equations, but 
that are useful in the discussion that follows are for the enthalpy h * 
e + P/p * CpT and for the difference Cp - c y * R/m and ratio Cp/c v * y 
of the specific heats at constant pressure and volume. 

To these equations must be added appropriate boundary conditions to 
represent properties of the incident solar wind and the planetary body. 
These include , T*, v„, and B* in the incident solar wind and P 0 , T 0 , 
v 0 * § 0 , $ 0 , and ion composition at some magnetoionosphere reference 
level between the magnetoionopause and the planetary surface. 

The difficulties presented in the solution of these equations are 
so great that neither analytic nor numerical solutions have yet been 
obtained for realistic three-dimensional boundary conditions 
representative of a planetary application. These difficulties arise 
from se/eral sources. First of all, the equations are nonlinear and the 
dissipative terms that contain the highest order derivatives in each of 
the differential equations are all multiplied by small coefficients. 
The shock waves and boundary and current layers that appear in the 
course of the solution are thus finite but thin regions of rapid 
variation at locations that cannot be specified in advance, but must be 
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found as part of the solution. Their adequate treatment requires the 
use of either fine spatial grids over extended regions or fine adaptive 

i 

grids In the vicinity of the surface; or the acceptance of poor 
resolution of the solution in these regions of rapid variations at best, 
and widespread contamination of the solution at worst. Use of fine 
spatial grids requires small time steps in unsteady simulations, or very 
many iterations tc converge a steady flow simulation. All of these 
factors contribute to the substantial computational requirements arising 
from the sheer magnitude of the problem. Eight quantities p, T, y, and 
B must be calculated at every grid point. For reasonable resolution of 
the total 3-D interaction region, there would be a minimum'of say 100 x 
100 x 100 * 10 6 grid points not allowing for extra grids in the vicinity 
of shock waves, boundary layers, etc*, and hundreds to thousands of 
iterations or time steps required to converge to a steady state solution 
from an assumed set of boundary and initial conditions. Obviously, the 
calculation of solutions will remain costly for some time, and we 
believe it will be many years before solutions of these equations will 
become inexpensive enough for extensive routine use in the analysis and 
interpretation of space observations. 

Howeyer, advances in computational capabilities are still 
increasing at a rapid pace, approximately doubling every 7 years 
according to Chapman (1979), due almost equally to advances in computers 
and in numerical methods. We may enquire at this point, therefore, as 
to what might be anticipated from the solutions of these equations if 
they could be attained. First, we believe these solutions would indeed 
provide good predictions of the bow wave location and the large-scale 
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features of the associated flow. Shock waves and magnetolonopause 

regions would have finite thickness, and magnetic merging and 

» 

reconnection would occur at the magnetolonopause. There would also be 

extensive heating and rarefaction at these boundary regions, 

♦ 

particularly those of the ionopause type, with the possible formation of 
magnetic barrier regions of high magnetic field separating the solar 
wind and ionospheric plasmas. 

We believe that these predictions would be correct qualitatively, 

but the quantitative predictions would be unreliable because of 

'inadequancies of the scalar transport coefficients. As noted previously 

(Spreiter, 1976), the expressions developed for nonmagnetized plasmas 

are highly inappropriate. For example, the widely used expression P - 

10” 1 6 T 5 / 2 « 0.469 y t p g cm” 1 s" 1 for the viscosity of 

P 

fully-ionized hydrogen having a representative value of 22 for the 
Coulomb logarithm, and in which T is the temperature in K, v t is the 
thermal velocity of the protons, p is the density, and is the 
effective mean free path for cumulative deflection of 90° by Coulomb 
interactions, leads to a Reynolds number R e = pvD/y of only 0.002 when 
the number density of protons Np = 10 cm" 3 , T = 10 5 K, and the radius of 
the Earth 6.37 x 10 8 cm is taken for 0. Such a value for R e is not at 
all indicative of an aerodynamic-like solar-wind flow past the Earth; 
but is more typical of a small ball sinking through tarJ Use of such a 
value for y would lead to the prediction of enormously thick boundary 
layers and shock waves, completely different from those observed. 
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There is no dilemma, however, since the particles were assumed In 

the derivation to travel in straight lines between collisions which, 

» 

even in the more conservative sense of cumulative small Coulomb 
deflections, are indicated to be separated by mean distances of the 
order of half an astronomical unit when the above-stated conditions are 
applied. This is obviously grossly inappropriate for planetary 
applications which involve phenomena of much smaller scale. 

Part of the answer to this apparent deficiency is provided by the 
fact that the presence of a magnetic field in the solar wind prevents 
the particles from traveling in straight lines between collisions, and 
causes them to spiral along the moving magnetic field lines. This . 
reduces the transport transverse to the field lines approximately as the 
square of the ratio of the gyroradius of the protons to the distance^. 

If a representative value of 1.4 x 10 " 6 is used for this ratio, 
corresponding to a magnetic field of 5 nT, the Reynolds number in the 
'example cited above would increase to 9 x 10 8 . Such a value is of the 
order of that encountered in ordinary aerodynamics, and is consistent 
with the presence of relatively thin shock waves and magnetoionopause 
surfaces and the generally good agreement between observations and the 
results of dissipationless fluid theories. 

However, all is not that simple. The magnetic field does not 
reduce the transport coefficients equally in all directions; in fact, 
it does not reduce the values for transport parallel to the field lines 
at all. The dissipative part of the associated model is thus highly 
anisotropic. The direction of the anisotropy, moreover, depends on the 
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magnetic field, and hence also the flow, and cannot be specified in 

advance. Since there is at present virtually no theoretical development 

* 

for the behavior of such a fluid for any application, the space 
scientist studying these features of the flow in terms of an anisotropic 
dissipative fluid is faced with the task of achieving major theoretical 
adyances or, as is more often the case, being satisfied with hopefully 

describing what he thinks will happen in qualitative terms based on 

analogy with the known behavior of isotropic fluids. In view of *he 

extreme anisotropy of the solar-wind plasma and the fact that the 
Coulomb deflection times upon which the analysis is based are much 
longer than the times required for solar-wind particles to traverse the 
significant part of a planetary flow field, it is evident that - 

considerable caution should be exercised in relying on such 

descriptions. It is, moreover, quite possible that this particle 
approach to transport coefficients may be only marginally relevant to 
the actual dissipative processes; and that the significant mechanisms 
are associated more with larger scale eddying and turbulent features of 
the flow as in many analogous situations with boundary layers in 
ordinary fluid mechanics. 

Finally, we may anticipate that the presence of time derivatives 

and nonlinear convective terms in equations (2) through (.9) almost 
certainly leads to fluctuations (turbulence?) in the solutions, 

particularly downstream of shock waves and with increasing distance from 
the nose along the magnetoionopause. These may, in fact, be correct 
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predictions, but their presence and snail scale introduces additional 

complications in the numerical solution that probably exceed by a 

* 

substantial margin present computational capabilities to resolve , 

Actual solution of the dissipative magnetohydrodynamic equations 
for planetary flows is still in its infancy (see Walker (1983) for a 
recent review). At the present time, about the only solution that has 
been worked out by using physically realistic, as opposed to numerical, 
dissipation is that of Brackbill (1982) illustrated in Figure 2. This 
is for unsteady two-dimensional, as opposed to three-dimensional, flow 
about a magnetic dipole, and is intended to demonstrate the sequence erf 
events which follow the southward turning of the Interplanetary magnetic 
field. The time units are Earth radii divided by the solar wind speed. 
As the southward field reaches the dayside magnetopause, reconnection 
begins at t = 35. Reconnection in the tail occurs when the rotational 
discontinuity reaches the nightside at t = 45. A magnetic bubble forms 
in the 'tail (n) which propagates down the tail. Earthward of the 
reconnection region, tail field lines first are stretched out (n^) and 
then snap back to a more dipolar configuration as tail reconnection 
begins. To verify that it is physical resistivity and not numerical 
resistivity that is responsible for the tail reconnection, Brackbill 
also ran his code with the electrical resistivity set to zero. In this 
case, reconnection occurs on the dayside because of numerical 
resistivity, but none occurs on the nightside. It should be recognized 
that these and other unsteady two-dimensional models are fundamentally 
deficient in that there can be no return flow around the Earth back to 
the dayside. Therefore magnetic flux must be returned to the dayside 
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artificially. In addition, questions of stability may have very 

different results in two and three dimensions. To illustrate, nobody 

% 

would marvel at a tight-rope walker in a two-dimensional world! On the 
other hand, the appeal of two-dimensional models is great in numerical 
simulations because of the substantial reduction in the number of 
numerical values involved. For the example cited above with 100 grid 
cells in each direction the reduction is from 8 x 10 6 to 6 x 10 4 values 
required to define conditions at any step in the calculation. 


’ DISSIPATIONLESS MAGNETOHYDRODYNAMIC THEORY FOR 
PLANETARY BOW WAVE FLOWS 

Because of the difficulties noted above, most analyses and 
discussions of planetary bow wave flows have been based on the 
dissipationless magnetohydrodynamic theory. This avoids altogether the 
uncertainties associated with the proper form of the dissipative terms', 
and all the mathematical difficulties associated with the thin, but 
finite, regions of rapid variations representing the bow wave, 
magnetoionopause, etc. The equations defining this theory are just 
equations (1) through (9) with all the terms involving u, S, k, and a 
set equal to zero. 

Let us examine the consequences of omitting the dissipative terms. 
There is no diffusion of vorticity, heat, or magnetic field*, there is 
no magnetic merging nor recombination, and magnetically separated 
regions can occur. The magnetic flux through a circuit that moves with 
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the fluid Is constant; giving rise to the "frozen field" concept. 
Entropy, and hence P/p y , is constant for a fluid element, and thus 
constant along a streamline for a steady flow, except when it crosses a 
shock wave. Since the amount of entropy increase varies along the bow 
wave, the entropy differs from one streamline to the next. Such a flow 
is often called isentropic, but not homentropic. This entropy condition 
reduces the number of dependent variables from 8 to 7 at each grid 
point. Shock wayes, boundary layers, and current sheets used to 
represent the bow wave, magnetoionopause, and other thin layers become 
mathematical discontinuity surfaces of zero thickness across which the 
plasma properties change discontinuously in definite ways governed by 
the conservation laws. This requires use in the numerical solution of 
either shock fitting techniques, probably with adaptive grids; or shock 
capturing methods with artificial viscosity, thermal , and electrical 
conductivities, either explicit or numerical, which thicken the 
discontinuity surfaces over several grid spacings. It is anticipated 
that most of the differences between these solutions and the corresponding 
dissipative ones are confined to the regions near these discontinuities; 
and that there is probably little change in most of the large-scale 
features of the predicted flow and magnetic field properties. Finally, 
we remark that the dissipationless equations are still difficult and 
costly to solve; and that very few solutions have actually been 
determined for realistic three dimensional flows at this time. 

There are, however, a number of definite and exact solutions of 
these equations for simpler boundary and initial conditions that are 
fundamental to both the qualitative interpretation of planetary bow wave 
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flows and the establishment of quantitative predictive models based on 
the dissipationless magnetohydrodynamic theory. Foremost is the normal 

t 

propagation speed C n of a plane waye. In gasdynamics, this is equal to 

— l/2 

the speed of sound C s and the resulting expression C $ = (vRT/m) ' 
follows directly from the dissipationless forms of equations (1) through 
(9) with the magnetic field B set to zero. This speed is the same in 
all directions, and can be represented in a planar cross section by the 
circle shown in the upner left of Figure 3 in which C s for a plane sound 
wave is illustrated as a function of angle 0 between the wave normal and 
an arbitrary coordinate axis. Sound wayes depend on the compressibility 
of the medium, and do not occur in an -idealized incompressible fluid. 
In magnetohydrodynamics, however, wave propagation can occur even in an 
incompressible conducting fluid. Such waves are the Alfven waves, and 
the diagram in the upper right of Figure 3 illustrates how the normal 
propagation speed of such a plane wave varies with direction such 
that C n = C A cos 0, where C A = ( B 2 / (4 ttp) ) 1 / 2 is the Alfven speed and 0 
is the angle between the wave normal and the magnetic field. The two’ 
lobes in this figure are thus circles osculating at the origin with 
their line of centers along the direction of the magnetic field. Plane 
Alfven waves can propagate in any direction, except exactly normal to 
the magnetic field, but it can be deduced from the properties 
illustrated that the disturbance created by a point source will 
propagate only along the direction of the magnetic field and with the 
Alfven speed C A . Values for C s and C A for fully ionized hydrogen plasma 
are shown in Figure 3, from which it can be concluded that both C<j and 
C A have values of about 25 to 50 km s" 1 for conditions representative of 
those in the solar wind at the orbit of Earth. 
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The propagation speeds of a magnetohydrodynamic plane wave In a 
dissipationless perfect gas are substantially more complex than either 
of the two limiting cases just discussed. There are, first of all, 
three distinct wave modes each with propagation speeds varying with 
angle 0, Of these, the rotational wave is simplest to discuss. It has 
exactly the same propagation speed and properties as the Alfven wave for 
an incompressible fluid. The other two plane wave modes, with 

c n = ± j{C 2 + Cjj ± [(C| + C 2 ) 2 - 4C|C^ cos 2 0] J / 2 } 1 / 2 (10) 

are called the fast and slow waves because their propagation speeds are 
greater or less than that of the rotational wave. The fast wave, which . 
has a plus sign before the inner square root, is of greater interest 
because it determines the region of influence and dependence in the 

solution, and in the flow it represents. Its normal velocity C n ^. varies 
from a maximum of (C s 2 + C^ 2 ) 1 / 2 in the direction perpendicular to the 
magnetic field to a minimum equal to the larger of C s or in the 

direction of the magnetic field. From the two bottom diagrams of Figure 
3, it may be concluded that the gasdynamic sound speed C$ approximates 
reasonably the magnetohydrodynamic fast wave speed for C s > C^, but 
poorly for C s < C^. 

For steady flows with velocity in excess of the fast wave speed, 
the Mach cones from an infinitesimal point disturbance are of prime 
interest. Their traces in the plane of the v and B vectors may be 

constructed either from the curves at the bottom of Figure 3, which are 
called Friedrichs I diagrams, by constructing the envelope curves called 
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the Friedrichs II diagrams and drawing tangents from the end of the 
velocity vector terminating at the origin (see e.g., Jeffrey (1966)); 

i 

or directly from the Friedrichs I diagrams by the procedure illustrated 
in Figure 4 (Spreiter et al., 1966b), In the latter, a circle is drawn 
with the velocity vector as a diameter, and the Mach cone is determined 
by the straight line from the end of the velocity vector to the 
intersection of the circle and a C n curve of the Friedrichs I diagram. 
There are thus, in general, three Mach cones, each associated with one 
of the three kinds of waves; and they are appropriately called the 
fast, slow, and rotational Mach cones. They correspond to the single 
Mach cone of gasdynamics shown in the upper right of Figure 4, which 
makes an angle w = arcsin (Mg" 1 ) with the velocity vector. The 
maginetohydrodynamic Mach angles depend not only on M s and M^, but on the 
angle 0. To illustrate the latter effect, four limiting cases are 
presented in the lower parts of Figure 4. For a given Mg and M^, the 
maximum and minimum Mach cone angles are given by co max = arcsin (Mg“ 2 + 
M a " 2 ) 1//z and the larger of u> m .j n = arcsin (Mg -1 ) or = arcsin (M^” 1 )* 
For B parallel to v , to,, = arcsin {.[(Mg 2 * M A 2 - l)/(Mg 2 + M^ 2 ) ] 1 / 2 } , and 
for B perpendicular to v, o»x = arcsin{((Mg 2 + M A 2 + 1 + [(Mg 2 + M A + l) 2 - 
4M$ 2 M^ 2 ] 1//2 }/(2Mg 2 M^ 2 )) l / 2 }. Lack of appreciation of these differences is 
widespread in the space research literature. 

The wave speed and Mach cone diagrams of Figures 3 and 4 are for 
infinitesimal disturbances, and are therefore not appropriate for 
dealing with the finite disturbances ass^ated with the planetary bow 
wave and any other shock waves that may deveiop in the flow. These are 
governed by the magnetohydrodynamic discontinuity relations, which 
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follow directly from the conservation laws and auxiliary relations of 

equations (1) through (9). With the unit vectors normal and tangential 

» 

to the discontinuity surface defined by n and t, the differences between 

any quantity Q downstream and upstream of the discontinuity designated 

by £Q] - Q 2 - q,, and v n - v nflQw - v nd1sc , which for steady flow is v„ 

= v n „. , the conservation laws for mass, momentum, and energy are 

''flow 
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These must be solved together with the Faraday magnetic jump equations 

<!t 5 n ‘ W * 0 ' “V ■ 0 < 14 > 

and the auxiliary condition for entropy C s] >. 0. In gasdynamics, these 
conditions with B set to zero completely determine the shock jump 
relations; but in magnetohydrodynamics an additional condition is 
required. It is the evolutionary condition (see, for example, Jeffrey 
and Taniuti (1964)), which states in essence that there must be a 
mechanism by which a succession of infinitesimal waves can evolve into a 
finite discontinuity. This distinction does not have to be made in 
gasdynamics because the entropy condition implies the evolutionary 
condition and conversely; but this is not so in magnetohydrodynamics. 
In the latter, the evolutionary condition implies the entropy condition, 
but the entropy condition does not imply the evolutionary condition; 
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and both must be satisfied. An example of a situation in which the 

entropy condition, but not the evolutionary condition is satisfied is 

* 

the case In which the shock normal Is parallel to the magnetic field and 
M s and M A are within certain ranges. Theory Indicates that a finite 
magnetohydrodynamic shock wave cannot eyolve from an accumulation of 
infinitesimal waves; and observational evidence Indicates that a well 
defined shock wave may indeed fail to develop under such conditions. 


Solutions of equations (11) through (14) show that fiye types of 
discontinuity surfaces may occur in dissipationless magnetohydrodynamic 
flows. They may b’e classified as two boundary types for which ? n = 


v nri - Vi. =0, namely 
n flow n disc 


tangential discontinuities 

B n = 0, [v t ] / 0, [B t ] t 0, [p] f 0, [P + B 2 /8tt] = 0 


(15) 


contact discontinuities 

B n f 0, [y] - [B] - [P] = 0, [p] f 0 


and three shock wave types for which v n i 0, namely 
rotational or Alfven shock waves 

v n = 9 = ±B n /(^7Tp) 1 / 2 , [v t ] = [B 1 ./(4 ttp) 1 / 2 ] 

[p] ■ [P] = [v n ] = [v 2 ] = [B 2 ] = [B n J = 0 
fast shock waves 

pv n > pv . tPl > O , [PI > 0 , [B n J = 0 , [B 2 ] > 0 
slow shock waves 

pv < pv , [p] > 0 , [P] > 0 , [B 1 = 0 , [B 2 ] < 0 
n np 11 


(16) 


(17) 


(IS) 


( 19 ) 
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All types may be anticipated In planetary flow fields if there are 
sufficient fluctuations In the solar wind or planetary flow boundaries, 

* i 

They may also be Improperly induced by growing perturbations In a 
numerical solution of the dissipationless equations, unless controlled 
by artificial viscosity and conductivities. 

More explicitly, a planetary bow wave must normally be a fast shock 
wave because pv n of the incident solar wind exceeds pv n ^ by a 
considerable margin over the forward part of the bow wave; and the 
magnetoicnopause must normally be a tangential discontinuity in order to 
separate the solar wind and planetary magnetic fields and plasmas. The 
outermost Mach wave of the McOn must be a fast Mach wave, and the . 
boundary between the flowing solar wind plasma and the lunar downstream 
cavity might be either a tangential or a contact discontinuity. In 
addition, all of these types of discontinuities can occur in the solar 
wind as a result of interactions of plasmas with differing speed and 
properties. 

The most completely worked out solution for dissipationless 

magnetohydrodynamic flow past a planetary obstacle is that of Spreiter 
and Rizzi (1974) for the special case of an axisymmetric obstacle in a 
flow in which B is parallel to v in the incident solar wind (see also 
Spreiter et al., 1970) for a related application to the Moon). Under 
these conditions B = Apy everywhere; and the dissipationless 

magnetohydrodynamic equations can be transformed into those of 
gasdynamics by introduction of new variables 
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V* * v £ X-X 2 p/ ( 4tt) ] , p* « p[i-X s p/(4TT)]“ l , 

•V »s* 


P* * P + B 2 /(8tt) , S* * S , 

« 

h* * h + [X 2 p/(4tt)] [l-X 2 p/(4ir)] , h * e + P/p (20) 

in which 

X 2 P /(4tt) * M a “ 2 - (p/pj M^ 2 (21) 

except that the equation of state P * p RT/m is replaced by a rather 
complicated relation that corresponds to no real gas. These 
pseudogasdynamic equations can be solved using the methods deyeloped for 
ordinary gasdynamics t to produce the results displayed in Figure 5. They 
show, for a given axisymmetric obstacle shape and fixed Mr * 10 and y * 

joo 

5/3, the location of the be./ wave for selected M* between 2.5 and 20. 
The conclusi; to be gained from these results is that the 
magnetohydrodynamic bow wave locations are virtually independent of M^ 
for values greater than about 10, and are practically identical with 
those of gasdynamics in which all effects of the magnetic field are 
disregarded. They also show that for small M Am the bow wave flares out 
more at the flanks, but approaches the obstacle closer near the nose. 
From this, it is evident that attempts to represent t'ne combined effects 
of M Soo and M^ by an "equivalent value" for M Sco to be used in a 
gasdynamic model cannot succeed. In gasdynamics there is a monotonic 
relationship between Me and distance to the bow wave in all directions 
from the obstacle. If Me is diminished, the bow wave moves farther 

Ooo 

from the obstacle everywhere. There is no way that an equivalent Mach 
number can be found that will simultaneously move the flanks of a 
gasdynamics bow wave farther outward and the nose farther inward. 
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Solutions for non-aligned fields are much more difficult to work 
out than for the aligned case, and very few results are available at the 

i 

present time. Among these Is the solution of Wu, Walker, and Dawson 
(1981) for three-dimensional, steady, magnetohydrodynamlc flow past a 
dipole magnetic field to simulate the Earth, for which constant density 
and pressure contours are presented In Figure 6. Although there is no 
magnetic field in the interplanetary plasma in their model, the results 
display the broad outlines of the observed flow. The results indicate a 
moderately defined bow wave across which conditions change approximately 
in accordance with the Rankine-Hugoniot jump conditions, and a poorly 
defined magnetopause. The locations of the magnetopause and bow wave 
resemble those observed in space, although it should be noted that the 

Mach number used in the calculations, i.e., M s a 2.5, is substantially 

00 

less than that normally observed in the solar wind upstream of the 
Earth's bow wave. Perhaps the least satisfactory aspect of the theory 
is the thickness of the bow wave and the magnetopause. The model 

boundaries are 4-5 Ea. ih radii (3-4 grid spaces) thick, whereas the 

observed boundaries are a few hundred kilometers thick at most, and 
usually considerably thinner. The calculated thickness of these 
boundaries results from the numerical dissipation in the simulation, 
which erroneously spreads over several grid spacings what should be a 

mathematical discontinuity surface in the dissipationless theory on 

which the model is based. It may be anticipated on the basis of 
continuing rapid improvements of computers and numerical algorithms that 
finer grids will be used in future models, thereby providing better 
resolution of these regions of rapid variation. As long as numerical 
dissipation is used in the calculations, however, caution should be 
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exercised In Interpreting the results, particularly those associated 

with real dissipative processes that may be expected to occur. Although 

* 

f he prediction of such features may tantallzingly resemble those 
actually believed to occur, it should be remembered that the associated 
quantitative values lack proper physical basis, and must be regarded 
more as suggestive than quantitative predictions. 


THE GASDYNAMIC APPROXIMATION OF MAGNETOHYDRODYNAMIC 

PLANETARY FLOWS 

In recognition of the great computational demands of full 

magnetohydrodynamic solutions, particularly as viewed in terms of the 
limited computing ability then available, an approximate theory of 
considerable utility was developed about two decades ago. Following on 
the success of purely gasdynamic calculations by Axford (1962), Kellogg 
(1962), and Spreiter and Jones (1363) of flow past an obstacle shape 
defined by the Beard (1960) and Spreiter and Briggs (1961) solutions of 
the classical Chapman-Ferraro theory of interaction of a collisionless 
stream of protons and electrons with a magnetic dipole, a more 

consistent approximate theory based on magnetohydrodynamics was 
developed by Spreiter et al. (1966a, b; 1967; 1968; 1969; and 

1970a, b ,c) , A1 ksne et al . (1967, 1970), and Dryer et al . (1966, 1967). In 
this theory, the magnetohydrodynamic tangential discontinuity jump 
conditions are retained for the boundary conditions at the magnetopause, 
and the bow wave is identified as a fast magnetohydrodynamic shock wave, 
which would tend to resemble a gasdynamic shock wave for large and 
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for M* > M s . The key approximation Is that the magnetic terms may be 

disregarded In the momentum and energy conservation equations on the 

* 

basis of their smallness for large M A . The resulting equations are thus 

00 

Perfect gas 

P = pRT/m (22) 


Conservation laws 
Mass 


It + 3x^ (pv k> = 0 

Momentum 


(23) 


3(pv i ) J . 3 , , n» , 9 i 9 k g*. . n 

TE + ITT ( p v,v k + P5 ik + to - ‘ sfe* - 0 


(24) 


Energy 


Pk (^— + pe + p$) + div[pv (Xj- + e + ~ + $)]= 0 


(25) 


and the Faraday Electricity and Magnetism Equations 

= curl (v x B) , di'v B = 0 (2G) 

The heat transfer equation degenerates to Ds/Dt = 0, which indicates 
that the entropy is constant along a streamline for dissipationless 
flow; but this relation does not apply across a shock wave because 
dissipation always occurs there in such a way that [ s] > 0 with the 

amount of entropy augmentation depending on the incident Mach number and 
the angle between the shock normal and the free-stream velocity vector. 
In most, but not all applications, the gravitational terms involving g 
and $ may be omitted. 
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With this approximation, the conservation equations become exactly 
those of gasdynamics. The magnetic field B therefore has no influence 
in the solution for the pressure p, the density P, the temperature T, or 
the velocity v. The latter quantities can be solved for using the 
methods of gasdynamics without further consideration of the magnetic 
field. The Faraday magnetic equations are the same as in the 
dissipationless magnetohydrodynamic theory, which implies that the 
magnetic field is "frozen" into the flow. However, the values for y In 
equation (21) are those of the gasdynamic solution. The resulting 
prediction for the magnetic field differs from that of an exact 

a 

magnetohydrodynamic solution by an unknown amount' dependent on the 
difference in the velocities in the two theories. This approximate 
separation of the magnetohydrodynamic problem into two parts, which can 
be solved in succession, leads to a significant reduction in computing 
requirements. The resulting component problems are still sufficiently 
complex that they have not been solved without further approximation. 
These difficulties are associated not only with the nonlinear and mixed 
el liptic-hyperbol ic type of the differential equations, but also with 
the discontinuity surfaces that represent the bow wave and the 
magnetoionopause whose locations are not known a priori but must be 
found as part of the solution. 

Problems associated with the unknown location of the 
magnetoionopause can be avoided to an acceptable degree over most that 
surface by use of an approximate relation for the pressure of the 
flowing solar wind on the magnetoionopause, i.e.. 
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P * K P„ v‘ cos 2 ip 


(27) 


This serves to decouple the calculation of the magnetoionopause shape 
from that of the surrounding flow. In this relation, K is a constant 
usually equated to unity and is the angle between the free-stream 
velocity and the exterior normal to the magnetoionopause with directions 
defined so that ^ s 0 at the magnetoionopause nose. With the proper 
selection of a value for K, this relation provides good results for the 
pressure over most of the magnetoionopause, but deterioration sets in as 
V approaches tt/ 2. For 'P >_ ff/2, equation (27) is totally inappropriate, 
and should not be used. 

Equation (27) has two distinct origins in planetary studies. On 
the one hand, there is the original usage in a long series of works 
starting with Chapman and Ferraro (1931) and particularly Ferraro 
(1952) in which the relation P = K m p n p v£ cos 2 ^ is derived by 
considering the interaction of a collisionless stream of protons of mass 
nip and number density N p and an equal number of electrons, all with 
uniform initial velocity y*,, with a dipole magnetic field oriented 
normal to the flow. Depending on the assumptions, K was usually found 
to be either 1 or 2, although 1/2 also appears occasionally in the 
literature of the time. That analysis was developed from a cold particle 
stream viewpoint to explain how a cavity that we would now call the mag- 
netosphere forms in a collisionless plasma stream or cloud of solar ori- 
gin that sweeps past the Earth and produces a geomagnetic storm. 
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The second line of argument that leads to equation (27) stems from 
the recognition that a highly supersonic solar wind exhibiting continuum 
fluid-like properties is present at all times; and that one of its 
manifestations is the bow shock wave that forms upstream of the 
magnetosphere, which acts much as an impermeable obstacle in the flow. 
From this point of view, equation (27) may be regarded as a close 
approximation for high Mach number flow to the Newtonian pressure 
relation for hypersonic flow P - - K p ro v£ cos z \p. This is usually, 
but not always, approximated in planetary applications by equation (27) 
in which P M is disregarded on the basis of its relative smallness as in- 
dicated by P = p v 2 /(yM 2 ). The constant K is determined from the Ray- 
leigh pitot formula for the stagnation pressure for supersonic flow 
decelerated to rest after passing through a normal shock wave. Under 
these conditions, K is found to be 



(y+1 )/ (y-1 ) I 


(28) 


For large Mach numbers, K = 0,881 for y = 5/3 and 0.844 for y = 2, the 
two most- commonly used values for the ratio of specific heats y in 
planetary flow studies. In using these pressure relations, it should be 
obseryed that equation (27) is not an exact relationship within any 
modern formulation of planetary flow theory. The original collisionless 
plasma beam theory is inappropriate because it disregards completely the 
deflection of the flow by the bow wave and in the intervening region; 
and the gasdynamic derivation is approximate because the Newtonian 
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pressure relation is a semi -empirical relation that Is not exact In any 
limiting case. That equation (27) does Indeed lead to good results 
compared with an essentially exact numerical gasdynamic solution for a 
number of representative solar wind conditions has been shown by 
Spreiter et al. (1966b, 1968). It seems to us equally plausible to 
consider for magnetohydrodynamic applications that Kp M ^ cos 2 ip should 
be -equated to the sum of the gas and magnetic pressures P+B 2 /8 tt or (P+ 
B 2 /87 t)-(P oo +B co 2 /8tt) rather than to P or (P-Pj alone as is customary. 


With the latter interpretation and disregarding P eo + B 00 2 /8iT as small, 
the tangential discontinuity boundary condition to be applied at a plan- 
etary magnetoionopause becomes 

Kp ro v 2 cos 2 i p = (P + B 2 /8tt) mip , B r = 0 (29) 

where subscript MI P refers to conditions just below the magnetoionopause. 
This enables the calculation of (a), the magnetoionopause location, and 
(b) conditions in the magnetoionosphere without further consideration of 
the solar wind plasma flow; and (c) reduces these problems for an 
idealized magnetic planet to those of the classical but physically 
outmoded Chapman-Ferraro theory (1931). 

We may summarize the specifications for such an idealized problem, 
and also its counterpart for a nonmagnetic planet, as follows. In 
either case, the normal component of the magnetic field is zero, and the 
solar wind pressure K p ro \£ cos 2 tp on the magnetoionopause is balanced by 
the sum of the magnetic and ionospheric gas pressures on the inside of 
the magnetoionopause in accordance with equation (29). Inside the 
magnetosphere of a magnetic planet such as the Earth, the intrinsic 
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planetary field dominates until it is terminated by the currents flowing 

in the magnetopause. Other current systems are present and important in 

% 

the calculation of magnetospheric processes; but, with the exception of 
the magnetotail current sheet which flows across the tail in 
approximately the geomagnetic equatorial plane, they are usually 
relatively unknown and unimportant, and disregarded in the calculation 
of the magnetopause shape and the surrounding flow. The field inside 
the magnetosphere is therefore just that of the intrinsic planetary 
field, usually approximated for the Earth by a magnetic dipole as 
distorted by the effects of the currents in the magnetoionopause and the 
tail plasma sheet. That field is calculated by solving 

div B = 0 , curl B = 4tj/c ( 33 ) 

subject to the boundary conditions 

Kp«, v £ cos ^ = ( B t/ 8tf ) M P » B n = 0 (31) 

at the magnetopause due to a planetary magnetic dipole field 

0^ = -(Mp/r 3 )( 0 _ sin 9 + r 2 cos 0 _) (32) 

at the origin. In this equation, M p is the planetary magnetic moment 
and r and 9. are geomagnetic radius and colatitude. In most applications 
j = 0 except at the magnetopause and in the tail current sheet. With 
neglect of the latter, the problem is exactly the Chapman-Ferraro 
problem. 
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For a nonmagnetic planet with a sufficiently dense ionosphere to 

stop the solar wind, the corresponding problem is specified by assuming 

» 

the ionospheric pressure above some reference level P Q can be 

approximated by 

P = P Q f(r,e) , B = 0 ( 33 ) 

within the ionosphere; and 

Kp M V* cos2] P = p ip (34) 

at the ionopause. In the original presentation of this model (Spreiter 
et al ., 1970a, b), the ionospheric pressure was assumed to be simply P « P 0 
expI-(r-r 0 )/H 1 where H = RT/mg * const, i.e., a constant scale height 
ionosphere was assumed. Subsequently, Spreiter and Stahara (1980a, b) and 
Stahara et al. (1980), included the inverse square variation of the 
gravitational acceleration g with r, while retaining the simple 
assumption of a constant temperature ionosphere. With improved 

knowledge, as has since been acquired for Venus, the only known clearly 
nonmagnetic planet, refinements in the specification of f(r,B) can 
obviously be made and incorporated into the calculations. We believe 
that little change in the ionopause shape will result, however, because 
the rate at which P diminishes with r is so great that only small 
changes in the ionopause shape will result. In other words, a 
nonmagnetic planet with a substantial ionosphere, such as Venus, 
presents a rather firm obstacle to the solar wind flow compared with a 
strongly magnetic planet as the Earth or Jupiter for which considerable 
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change in the size of the magnetosphere occurs with variations in the 
solar wind momentum flux p vi. 

' oo 00 

» 

Figure 7 illustrates the resulting shape calculated for the 

magnetopause associated with a magnetic dipole of magnitude and range of 
orientations in the solar wind resembling that of the Earth. These 
results are for the first-order approximation of Beard (1960) as solved 
by Spreiter and Briggs (1962) and Briggs and Spreiter (1963). A number 
of higher-order solutions have been worked out subsequently, but the 
results are -virtually indistinguishable on the scale shown and the 
differences are not particularly relevant to the present discussion. 
(See Wu and Cole, 1982, and Wu, 1983, for recent commentary.) These 
results, and all others of the kind, have been based on the 
Chapman-Ferraro formulation of the interaction. However, they may 
equally be considered to be solutions of the continuum gasdynamic or 
even magnetohydrodynamic models to the extent that the solar wind 
pressure on the magnetopause can be approximated by the extended 
Newtonian relation of equation (29). From the latter point of view, 
inaccuracies in the calculated magnetopause shape may be anticipated 
along the magnetosphere tail and near the neutral points because the 
Newtonian pressure relation becomes increasingly inaccurate as ip 
approaches tt/ 2. The neutral points are of particular interest because 
they represent points in the idealized theory where the magnetospheric 
magnetic field vanishes and particle entry into the magnetosphere may 
occur. There is considerable evidence that such entry does indeed 
occur, but any quantitative prediction of conditions in the vicinity of 
the neutral points that uses the Newtonian pressure relation in the 
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boundary condition as given in equation (31) or (33) must be regarded as 
unreliable. The reason is that the condition B * 0 at a neutral point 
leads through the boundary condition to cos 2 t|; * 0, or * tt/2; but that 
is precisely the condition for which the Newtonian approximation clearly 
fails. In the Chapman-Ferraro plasma beam approach, a further 
deficiency is that effects of double impacts, as would occur near a 
neutral point when a particle grazes the magnetopause just forward of 
the neutral point and hits the magnetopause again just aft of the 
neutral point, are not considered. In the continuum gasdynamic or 
magnetohydrodynamic models, the presence of an indented region results 
in the formation of an embedded shock attached to the surface (Walters, 
1966); but Spreiter and Summers (1967) pointed out that such a shock 
could not develop adjacent to the magnetopause. The reason is that the 
pressure jump across the shock wave cannot be matched with a 
corresponding jump in the magnetic pressure in the magnetosphere. The 
resolution they proposed was that the flow would separate somewhat ahead 
of the neutral point and then reattach downstream of the neutral point, 
and that the constant pressure along the intervening free streamline 
surface would cause a cusp-shaped plasma filled region to form. Recent 
magnetohydrodynaraic numerical solutions of Wu (.1983) support this 
proposal. With the cusped geometry, it is no longer required for the 
magnetic field to vanish at the netural point, but solar wind particles 
can still enter the magnetosphere near the tip of the cusp. 
Considerable interest has been attached to observations in these polar 
cusp regions in recent years, and also to the conditions in the mantle 
region which extends downstream from the cusped regions along the 
magnetopause. 
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From the present point of view of planetary bow waves and their 
associated flow fields, perhaps the main point is that there exists a 

i 

well worked out body of theory for the calculation of the magnetopause 
shape and location, and that it has been found in a great many studies 
to provide a reasonably good prediction of conditions actually observed 
in space. This is true for both the Earth and for other planets for 
which data are available, provided allowance is made for the differing 
nature of the planetary obstacles. For example, the ionopause shape for 
the nonmagnetic planet Venus is axisymmetric in the theory and nearly so 
in observations; and there are of course no neutral point phenomena. 
Jupiter's magnetosphere appears to be flattened somewhat toward the 
equatorial plane, which is understandable in terms of the effects of its 
extremely large trapped particle regions and high spin rate. 

To proceed with the solution of our problem, the next step is to 
calculate the supersonic gasdynamic flow around the magnetopause, which 
can at this stage be considered to be similar to an impermeable obstacle 
of known shape. For a nonmagnetic planet, this obstacle is axisymmetric 
about the line through the planet center directed along the free-stream 
velocity direction. For a magnetic planet, as the Earth, however, the 
obstacle shape is nonaxisymmetric and a decision has to be made whether 
or not to include the effects of the moderate departure from 
ax i symmetry. When the original calculations were made over twenty years 
ago, there was no choice. Only the theory for axisymmetric flow had 
been developed, and even that was very new at the time. The computers 
and methods then available were totally inadequate for dealing with a 
fully three-dimensional supersonic flow past a blunt body. As a result, 
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a further simplification was introduced in that the magnetopause was 

approximated by an axlsymmetrlc shape. This was usually, but not 

» 

always, taken to be the body formed by rotating the equatorial trace of 
the magnetopause surface. This choice was motivated partly by that fact 
that most space observations for the Earth were made near the equatorial 
plane. Today, adyances in computers and numerical methods, and recent 
motivation provided by the needs of the nonaxi symmetric space shuttle as 
opposed to the axisymmetrlc vehicles of the earlier era, have made 
available methods for calculating supersonic gasdynamlc flows about 
nonaxi symmetric blunt objects. However, the cost in computing time 
rises substantially for these 3-D flows, and the errors from assuming an 
axi symmetric magnetopause shape are anticipated to be no larger than 
those from a number of other sources in the entire procedure. 

Even with the assumption of an axisymmetric obstacle, the 
determination of the gasdynamic flow properties is the most difficult 
and time-consuming portion of the numerical solution of the entire 
planetary flow and magnetic field problem. As sketched in Figure 8, two 
different methods are used in the calculations; one for the nose 
region, where both subsonic and supersonic flows occur, and another 
downstream of that region, where a more computationally economical 
procedure can be employed since the flow is supersonic everywhere. In 
the earlier analyses, an inverse iteration method was used for the nose 
region, and the method of characteristics was employed for the remaining 
supersonic region. Both of these procedures are inferior from a 
computational point of view to more' recently developed methods and were 
replaced in the modernization of these solutions described by Spreiter 
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and Stahara (1980a,b) and most comprehensively by Stahara et al. (1980). 

The nose -region is treated using a new ax i symmetric implicit unsteady 

» 

Euler equation solver specifically developed for the present 

application. That procedure determines the solution in the nose region 
by an unsteady asymptotic time-marching procedure which advances the 
solution forward in time until the steady state is attained. The 
remainder of the flow field is determined by a shock capturing marching 
procedure which spatially advances the solution downstream as far as 
required by solving the steady Euler equations. Inasmuch as these 
methods are complex, lengthy, and fully described in the above 
references, we turn directly to a discussion of the results. 

Figures 9 and 10 show typical results for the Earth for 

representative solar wind conditions of y a 5/3 and M s * 8. In addition 

00 

to the bow wave and the magnetopause, Figure 9 shows the sonic line 
which delineates the boundary between subsonic and supersonic flow, the 
streamlines which indicate the local direction of the flow, and the Mach 
or characteristic lines which indicate the zones of influence and 
dependence of an infinitesimal disturbance in the flow. Figure 10 shows 
contour lines of constant density, velocity, and temperature. These are 
sharply defined solutions, with much finer resolution than the 
magnetohydrodynamic solutions presently available, and can be determined 
in one to two minutes on a large modern computer as opposed to several 
hours for a much less well resolved magnetohydrodynamic solution. An 
extensive catalog of results for various interplanetary conditions is 
given in Stahara et al. (1980). 
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With the flow field determined by the gasdynamie calculations, the 

magnetic field B may be determined by integrating equations (14) and 

* 

(26); or alternatively the following equations derived from them 

wl h S • “5 - 0 - -jhs • V l < 35 > 

in which S is an arbitrary surface moving with the fluid and D/Dt * B/3t 
+ (v • v) is the material or substantial derivative. As noted 

*■* »v 

previously, these equations are commonly Interpreted as indicating the 
field lines move with the fluid. For the steady state in which 3/H * 0 
and v n 3 v n , these equations lead to a straightforward calculation in . 
which the vector distance from each point on an arbitrarily selected 
field line to its corresponding point on an adjacent field line in the 
downstream direction is determined by numerically integrating /v dt *As 
over a fixed time interval At. Once the magnetic field lines are 
determined, the magnetic field at any point may be calculated from the 
relation 

B/B* ■ (p/pJ(AJ/AJlJ (.36) 

where M is the vector length of a small element of a flux tube. In 

•V 

closing this discussion, one should observe that, in general, both 
streamlines and magnetic field lines are three-dimensional curves. For 
axisymmetric flows, the streamlines are planar curves, but the field 
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lines are three-dimensional curves except fov the special case of 
aligned field flow in which the B and v vectors are locally parallel 
throughout the flow. 

The foregoing procedure is valid generally, but for axisymmetric 
flows It is much simpler and completely equivalent to compute the 
magnetic field by a decomposition theorem developed by Alksne and 
Webster (19/0) whereby the magnetic field at P is given by 



can be calculated directly from the gasdynamic solution. (B p /B ro ) x must 

still be calculated using the As - Ail method, but only for the plane of 

magnetic symmetry in which 0 and for which B p is two-dimensional. 

In this decomposition, the subscripts ll , ± , and n refer to 

contributions associated with the component B^ of B w parallel to v OT ; 

the component B OTi perpendicular to v w in the plan i that contains the 

point P, the center of the planet, and the vector v M ; and the component 

normal to the latter plane, where e„ is a unit vector in the latter 
n n 

direction. The quantity r p in equation (38) is the radial cylindrical 
coordinate of the streamline through P. Figure 11 is included to 
facilitate understanding of the quantities involved in this analysis, a 
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key point of which is the orientation of the coordinate system so that 
the point P at which the magnetic field is to be evaluated is in the xi. 

i 

plane. 

Representative results for the magnetic field direction and 
intensity are shown in Figure 12 , The upper part of the figure presents 
the magnetic field lines and contours of equal field magnitude in the 
plane of magnetic symmetry, for which the magnetic field lines are all 
planar calves. Out of this plane, the field lines are three-dimensional 
curves and two projections are necessary to convey the results. These 
magnetic field results vary greatly with changes in the interplanetary 
field direction, but are relatively insensitive to moderate variations 
in the other flow parameters. 

A property of the magnetic field calculated in this way that is 
sometimes cited as a deficiency of the theory, but which can actually be 
understood in terms of real processes as viewed in the limit of 
vanishing magnetic field and dissipation, is that the calculated field 
becomes very large near the magnetoionopause. This may be deduced 
directly from equation (36) for the general case or equations (37) and 
(38) for axisymmetric flows from the fact that BpP co /(B 00 P p ) approaches 
infinity as P approaches the magnetoionopause, because the length A£/ A£ ra 
of the element of a convecting magnetic flux tube directly incident upon the 
magnetoionopause nose enlarges indefinitely as it flows past the 

magnetoionopause. It should be noted that this is not just for the 
gasdynamic-convected- field model, but applies equally to the 
dissipationless magneto hydrodynamic theory since equation (36) is an 
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exact relation in that theory. The only difference between the 
gasdynamlc and magnetohydrodynamic procedures is the difference in the 
values for v in the calculation of As * / v dt. Such changes will 
inya idate the above statements only if v is zero in the 
magnetohydrodynamic solution over the entire magnetoionopause instead of 
just at the stagnation point and that is not the case. The theory thus 
indicates that either Bp/B^ is infinite or p p is zero. Both 
possibilities are meaningful in the limit of vanishing dissipation and 
interplanetary field as we explain below. 

If the boundary is of the pure magnetopause type with a vacuum 
magnetic field in the magnetosphere, the interpretation is that B D 

r 

increases to the finite magnetospheric level, which is independent of 
the interplanetary field B^; and p diminishes to zero as the point P 
moves across the magnetopause from the solar wind plasma to the 
magnetosphere. Now the theory does not provide values for Bp directly, 
but only for the ratio Bpp OT /(B^ p p ) and that strictly only in the limit 
of vanishing B^. The infinite value for this ratio is thus the theory's 
best effort to predict a finite value for B p that is independent of B^ 
as B ot approaches zero. In the dissipationless magnetohydrodynamic 
theory for which the limiting process of vanishing B^ would not be 
applied, the magnetic field would be expected to rise to the 
magnetospheric level in much the same way as the point P approaches the 
magnetopause. Meanwhile the density would drop in such a solution and 
become exactly zero at the idealized magnetopause; so that in that 
interpretation of the theory Bpp^/CBjDp) becomes infinite at the 
magnetopause, but for a different reason. This is, in fact, the 
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mechanism by which a magnetopause would form in dissipationless 

magnetohydrodynamic theory with distinct magnetic field configurations 

* 

on the two sides of the boundary for steady flows. 

If, on the other hand, the boundary is of a pure ionopause type for 
which there may be no ionospheric magnetic field, the prediction of 
infinite B p /(B^ p p ) may still be understood, although In a somewhat 
different way. As point P approaches the ionopause, the magnetic field 
would again rise, but this time the rise would be limited by the 
ionospheric pressure, which is independent of the interplanetary field. 
Thus the ratio of the finite ionopause magnetic field to the 
interplanetary magnetic field is infinite in the limit of vanishing 
interplanetary magnetic field. If the magnetic field actually rose to 
such a value as to balance the ionospheric pressure, the gas pressure 
would be driven toward zero. This would indicate the presence of 

magnetic barrier region of high field and little plasma separating the 
flowing solar wind and the ionospheric plasmas. From the higher point 
of view of dissipative theory, an ionopause is a shear layer in 
hypersonic flow in which intensive heating occurs. Since the pressure 
on the ionopause is determined primarily by the momentum balance and is 
nearly independent of the dissipative processes, a rise in temperature 
leads, through the perfect gas law P = p R T/m, to a fall in density in 
the shear layer. In this way, a barrier region of hot rarefied plasma 
would form between the solar wind and the ionospheric plasma even in the 
absence of all magnetic fields. From several points of view, therefore, 
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It appears plausible to expect these theories to predict enhancement of 

the Interplanetary magnetic field near a magnetolonopause and a 

» 

relatively evacuated magnetic barrier region at an lonopause. 

The above account provides an oyerview of the essential parts of 
the gasdynamic-convected-field model, and its relationship to the 

dissipative and nondissipative magnetohydrodynamic theories from which 
it is derived. Once the flow and magnetic field properties are 
determined in the -manner described, further extensions such as to the 
calculation of the electric field E = -(v x B)/c,. proton velocity 

<S# (V »w 

distributions, and many other quantities can be accomplished in a 
straightforward way with the introduction of the appropriate formulae. 
Because it is the most completely worked out method for calculating 
planetary bow waves and their associated flow fields and is gaining 
widespread use in the interpretation and analysis of data, we feel it 
appropriate in concluding this section to make a number of general 
remarks about the present state of this theory as implemented by the 
computer model of Stahara et al . (1900), and Spreiter and Stahara 

(1980a, b). 


First of all, it should be understood that the 
gasdynamic-convected-field model is essentially a magnetohydrodynamics 
theory in the limit-of vanishing dissipation and vanishing B OT , i.e., for 
M a +<*with > M s « Numerical solutions are obtained after the 

oo 

introduction of additional approximations, particularly the Newtonian 
pressure relation and an axisymmetric obstacle shape. The entire model 
is fully implemented with a documented computer code in NASA CR 3182 
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(Staharaet al. , 1980). An entire calculation requires about 120 seconds for 
a single case on a CDC 7600 or about 20 seconds on a Cray XMP, About 90 

i 

percent of this time is for the gasdynamic part of the calculation, with 
only about 10 percent being required for the magnetic field. The model 
is proving to be a useful tool in the interpretation of data, 
particularly in understanding effects of the three-dimensional nature of 
the interaction process. Several examples of this may be found in other 
articles in this conference proceedings. Extended applications 
currently underway include: 

• bow shock shape and position for all the terrestrial planets 
(Slavin et al ., 1981, 1983) 

• distant planetary Mach cone and bow shock shapes for Venus, 

Earth and Mars (Slavin et al., 1984a) 

• magnetospheric source of energetic particles upstream of 
Earth's bow shock (Luhmann et al., 1984b) 

• locations and asymmetries of magnetic field merging sites on 
the Earth's magnetopause (Luhmann et al., 1984a; Crooker et 
al., 1984a) 

• magnetic field draping ori the Earth's magnetopause (Crooker 
et al., 1984b) 

• intrinsic magnetic field of Mars (Russell, 1984) 
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* detailed plasma and field studies In the Venusian Ionosheath 
(Mlhalov et al., 1982) 

i 

• bow shock shape and position for Jupiter and Saturn (Slavin 
et al., 1984b) 


CONCLUSIONS 


To summarize the preceding discussion of magnetohydrodynamic and 
gasdynamic theories for planetary bow waves, we conclude with the 
following remarks.- 

Magnetohydrodynamic theory provides a good basis for modeling many 
aspects of planetary bow waves and their associated flow fields. 
Current development of numerical solutions for representative 
applications is in the early stages, and presently require extensive 
computer times. Even when these solutions are obtained, there remain a 
number of important phenomena that cannot be described by the theory. 
Examples include upstream waves or any other multi -component plasma pro- 
cess. Moreover, the dissipationless limit of the general MHD theory can- 
not account for any phenomena that depends on the finite thickness of 
shock waves, boundary layers, or current sheets. 

The aligned magnetohydrodynamic case in which the magnetic field is 
aligned with the velocity can be reduced to a pseudogasdynamics problem 
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and solved by methods of gasdynamics. The results display substantial 
effects of variations of the Alfven Mach number as well as the sonic 
Mach number, and are of such nature that they cannot be made Into a 
function of one parameter by any combination of the sonic and AlfvSnlc 
Mach numbers. 

The gasdynamic-convected-field theory is a limiting form of 
magnetohydrodynamic theory for large Alfven Mach number and for Alfven 
Mach numbers in excess of the sonic Mach number. With the aid of 
certain additional approximations including the Newtonian pressure 
relation for predetermining the magnetoionopause shape and the 
introduction of an axi symmetric magnetoionopause shape, thetheorycan 
be reduced to a form that is both amenable to efficient numerical 
solution and satisfactorily accurate for many purposes. This theory has 
been fully implemented in a readily available model that has been 
demonstrated to give useful results over a wide range of conditions. 
Several groups are now actively using the model in their work, 
developing further extensions and applications, and generally finding 
the results useful for both interpreting data and developing new 
insights. 

Finally, we note that the dissipationless, gasdynamic-convected- 
field theory contains certain properties that are sometimes interpreted 
as false and without meaning, but which actually can be understood as 
true indicators of real processes when viewed in the limit of vanishing 
dissipation and interplanetary magnetic field. A discussion of several 
examples is provided. 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6. 


Figure 7. 


Illustration of solar wind flows past major bodies w the 
solar system. 

Two-dimensional magnetohydrodynamic simulation of evolution 
of the geomagnetic field in a substream (Brackbi 11, 1932). 

Magnetohydrodynamic propagation speeds for plane waves. 

Magnetohydrodynamic Mach cone traces delineating regions of 
influence of a small disturbance at a fixed point in a 
steady flow. 

Calculated bow wave positions for various Alfvdn Mach 
numbers M^ in a steady aligned magnetohydrodynamic flow 

co 

with M_ -10 and y =5/3 (Spreiter and Rizzi, 1974). 

00 

Results of three-dimensional magnetohydrodynamic simulation 
of steady solar wind flow past the dipole magnetic field of 
the Earth in the absence of an interplanetary magnetic field 
(Wu et al ., 1981). 

Shape and size of the Earth's magnetopause according to the 
Chapman-Ferraro theory, or equivalently the Newtonian 
gasdynamic model. 


Figure 8. Outline of two alternative methods for calculating 

axisymmetric supersonic flow past a blunt-nosed body. 

* 

Figure 9. Calculated gasdynamic streamlines and Mach lines for steady 
supersonic solar wind flow past the Earth’s magnetosphere, 
M s =8, Y=5/3 . 

CO 

Figure 10. Calculated gasdynamic density, velocity, and temperature 

contours for steady supersonic solar wind flow past the 

Earth's magnetosphere, M- =8, Y=5/3. 

00 

Figure 11. Diagrams to illustrate the xin coordinates and magnetic 
field components used in Equations (37) and (38) to 
implement the Alksne-Webster, (1970 ) magnetic field 
decomposition theorem for axisymmetric flows. 

Figure 12. Magnetosheath magnetic field magnitude and direction for 
representative conditions for solar wind flow past the 
Earth, M«. =8, y =5/3 . 
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